Angular momentum conservation in heavy ion collisions at very high energy 
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The effects of angular momentum conservation in peripheral heavy ion collisions at very high 
qq". energy are investigated. It is shown that the initial angular momentum of the quark-gluon plasma 

should enhance the azimuthal anisotropy of particle spectra (elliptic flow) with respect to the usual 
picture where only the initial geometrical eccentricity of the nuclear overlap region is responsible 
, for the anisotropy. In hydrodynamical terms, the initial angular momentum entails a non trivial 

dependence of the initial longitudinal flow velocity on the transverse coordinates. This gives rise 
to a non-vanishing vorticity in the equations of motion which enhances the expansion rate of the 
supposedly created fluid compensating for the possible quenching effect of viscosity. A distinctive 
signature of the vorticity in the plasma is the generation of an average polarization of the emitted 
hadrons, for which we provide analytical expressions. These phenomena might be better observed at 
LHC, where the initial angular momentum density will be larger and where we envisage an increase 
of the elliptic flow coefficient V2 with respect to RHIC energies. 
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^ ' I. INTRODUCTION 

c, 

Nuclei colliding at ultrarelativistic energies have a large initial orbital angular momentum L$ if their impact pa- 
£C} . rameter is of order of some fm; in fact, for symmetric nuclei, Lq ~ Ayfs NN b/2 in natural units (H = 1). For Au-Au 
' collisions at RHIC energies v^atat = 200 GeV and Lq ~ 5 x 10 5 at an impact parameter b — 5 fm. The angular- 
momentum will be almost two order of magnitude larger in the forthcoming Pb-Pb collisions at LHC, at \/~s NN = 5.5 
TeV, with Lq ~ 1.4 x 10 7 . Due to the inhomogeneity of the colliding nuclei in the transverse plane, a significant 
fraction of L must be deposited in the interaction region, in other words should be transferred to the supposedly 
. ■ formed Quark-Gluon Plasma (QGP). Large values of the initial angular momentum of the plasma may give rise, as 
we will show, to significant observables effects. 

According to the to-date generally accepted description of the collision process, a locally equilibrated plasma is 
formed after a relatively short proper time (of the order of 1 fm/c) followed by a purely ideal-fluid hydrodynamical 
expansion. This kind of approach proved to be able to reproduce the large observed values of the elliptic flow in 
peripheral collisions, at a finite impact parameter, and the transverse momentum spectra of particles in the low px 
region |lj. Usually, in this kind of description, the Bjorken hydrodynamics scaling hypothesis is used either all along 
the evolution (2+1 hydro) Q or just at the initial proper time (3+1 hydro 0). In both cases, the initial longitudinal 
flow velocity only depends on z which amounts to make the initial angular momentum vanishing unless the energy 
density has an asymmetric dependence on the transverse coordinates [3j. But even if this was assumed, and the 
initial angular momentum was then non-vanishing, the dynamical evolution would be different from the case of a 
longitudinal flow velocity depending on transverse coordinates, as we will show later. 

In recent papers 0, H[ it has been found that amending the ideal fluid assumption with even a minimal viscosity 
strongly affects the elliptic flow. Particularly, Heinz and Song pointed out that in order to restore the agreement with 
a hydrodynamical description, one should enforce significant modifications of the initial conditions or the equation of 
state, such that the authors raise some doubts about the interpretation of RHIC results. In this paper, we want to 
show that including the initial angular momentum by a suitable modifications of the initial fluid velocity profile may 
cure the problem, or at least it may give a contribution in this direction. In fact, a finite angular momentum enhances 
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FIG. 1: Sketch of a peripheral heavy ion collision at very high energy in the longitudinal projection. The initial momentum 
distribution of the interaction region (right) should have a gradient along the axis x orthogonal to the collision axis z stemming 
from the different transverse densities of the colliding strips (left). 

the elliptic flow coefficient and broadens the transverse momentum spectra, exactly what is needed to counterbalance 
the quenching effect of viscosity. 

The most distinctive signature of an intrinsic angular momentum would be the polarization of the emitted hadrons. 
This argument has been put forward in refs. d, 0] where the authors take a QCD perturbative approach. Also, 
more recently, polarization has been related to the fluid vorticity Q , yet without developing an explicit mathematical 
relation. In this paper, we take advantage of a very recent study of the ideal relativistic spinning gas [9j and present a 
formula relating polarization to the angular velocity of an equilibrated, i.e. rigidly rotating, hydrodynamical system. 
We argue, on the basis of the locality principle, that such formula should hold for the most general fluid motion where 
the angular velocity is to be presumably replaced by an expression involving the local acceleration, hence the vorticity, 
of the fluid. 



In the usual picture of a peripheral heavy ion collision at ultrarelativistic energy the overlapping region of the two 
incoming nuclei gives rise to QGP whereas the non-overlapping fragments fly away almost unaffected. Thereby, only a 
fraction of the initial angular momentum Lq is left to the interaction region, while the largest part is carried away by 
the fragments (see fig.[T]). The angular momentum of the interaction region takes its origin from the inhomogeneity of 
the density profile in the transverse plane, the so-called thickness function. This is much clearly seen in a longitudinal 
projection: the colliding strips of nucleons have, in peripheral collisions, different number of nuclcons. While the 
central strips have the same weight, the strips above it will have a net momentum directed along the negative z axis 
and conversely the ones below it (see fig. [J). The net momentum density at each point (x,y) of the overlap region in 
the transverse plane (see fig. [5] for the axes definition) for symmetric (equal nuclei) collisions reads: 



II. ANGULAR MOMENTUM CONSERVATION IN HEAVY ION COLLISIONS 
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where T(x,y) is the thickness function, i.c 



the longitudinal integral of the nucleon density: 
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FIG. 2: Sketch of a peripheral heavy ion collision at very high energy in the transverse projection. The overlap almond-shaped 
region is marked in light grey and has an overall angular momentum directed along the symmetry axis y, orthogonal to the 
reaction plane. 



Only if the two colliding objects were homogeneous in the transverse plane, the angular momentum of the interaction 
region would be vanishing. Yet, the nuclei are not homogenous in the transverse plane; for instance, if they are 
assumed to be homogenous spheres in their rest frame, their thickness function T(x, y) would be proportional to 
2 , r being the distance from the centre of the nucleus and R its radius. In this case, eq. ([1]) would become: 
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From this momentum density, one gets an initial angular momentum J of the interaction region directed along the y 
axis: 



J = 2riQ dx dy x 



yjR 2 ~y 2 -{x- b/2) 2 - y/B? -y 2 -(x + b/2)' 



(3) 



In fig. [3]we show J for two colliding Gold nuclei at \fs NN = 200 GeV, in the two cases of hard spheres and Woods- 
Saxon distribution. For the former case, it is seen that the angular momentum attains a maximal value at an impact 
parameter of 2.5 fin and quickly drops thereafter. The maximal value of J is about 7.2 x 10 4 , i.e. 29% of the 
initial orbital angular momentum Lq of the colliding nuclei at that impact parameter. Therefore, J is very large and 
strongly dependent on the impact parameter b but this effect is usually ignored in the initial conditions assumed for 
hydrodynamical calculations as in the commonly used Bjorken model the longitudinal flow velocity only depends on 
z and it does not thus have any azimuthal anisotropy. 

This can be seen again from fig. [TJ since the net momentum of the colliding strips varies monotonically along x, 
either the proper energy density or the fluid four-velocity or both must have an asymmetric profile in x for the initial 
angular momentum to be conserved. If we take the reasonable assumption that the proper energy density cannot have 
such an asymmetric dependence on x because it can only depends on the density of nucleons at each point, the only 
remaining possibility is to admit that the initial longitudinal flow velocity is asymmetric in x from the very beginning, 
i.e. it is azimuthally anisotropic, in such a way that: 
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FIG. 3: Angular momentum J of the interaction region as a function of the impact parameter for Au-Au collisions at y/s NN = 
200 GeV. 

for a perfect fluid and if the initial flow transverse flow velocity is zero; in the above equation T is the stress-energy 
tensor, p is the proper energy density, p the pressure and 7 2 = (1— v )~ . The fact that v z is not azimuthally isotropicis 
implies, in general, a non-vanishing vorticity u> = (1/2) V x v for the fluid motion, and this may have remarkable 
consequences on the final particle spectra. It should be pointed out that some calculations || indeed introduce an 
x-asymmetric proper energy density function. Still, even if p was forced to have such an asymmetric x dependence 
in order to fulfill angular momentum conservation Q, the final velocity field would not be the same as when, more 
reasonably, v z is asymmetric in x. We will try to illustrate such effects with an oversimplified hydrodynamical scheme 
in the next section. 




III. HYDRODYNAMICAL SCHEME 



We are now going to set up a very simple hydrodynamical scheme to show that the azimuthal anisotropy of the 
longitudinal flow velocity required by the angular momentum conservation must enhance the elliptic flow. 

As has been mentioned, the requirement of an initial azimuthal anisotropy of the longitudinal flow velocity breaks 
the usual Bjorken scheme, where v z = z/t. As a first step, one would like to introduce a minimal change of the 
Bjorken scheme, which is not an easy task though. Thus, to describe the possible effects of an initial dependence of 
v z on x in the most transparent way, we will assume an oversimplified scheme assuming that the two colliding nuclei 
give rise to a complete thermalization within an infinitesimally thin slab Az at the time t = (see fig-H]). This scheme 
looks very similar to Landau hydrodynamical model, were not for the inclusion of an initial flow velocity v z (x) which 
ought to vanish in x — for an infinitely thin slab, for evident symmetry reasons. Such a picture of the collision 
should be more realistic at asymptotically large energies, where one expects thermalization to be extremely quick and 
nuclei are infinitely Lorentz-contracted along their collision axis. Furthermore, we will assume to deal with a perfect 
fluid and we will focus our attention on the transverse motion only. 

Firstly, we can write a relation for the initial momentum density 

/ ^ 2 1 dP 

(Pa+Poho v zO = ^-fafaj ( 5 ) 



where po,Po and v z q are the proper energy density, pressure and longitudinal flow velocity at the time t = 0, with 
7o = 1/(1 — v zo) because the initial transverse velocity is vanishing; dP/dxdy is given by eq. (fTJ) . The eq. © makes 



FIG. 4: Initial longitudinal velocity profile for the limiting case of sudden thermalization in the very thin overlap region of the 
colliding ultrarelativistic nuclei. 



it clear that v z q in this approach cannot be zero because of the initial unbalance in momentum and that, in general, 
depends on both x and y. The specific functional form of the proper energy density p affects the functional dependence 
of v zQ but it should not suppress its dependence on x because it is reasonable to assume that it has a symmetric 
dependence on x, unlike dP/dxdy, as pointed out at the end of Sect. 2. From a hydrodynamical point of view, the 
remarkable consequence of this is that the initial vorticity u) = (1/2) V x v is non- vanishing, unlike in the traditional 
Bjorken picture where v z only depends on z. Indeed: 

^ (t = 0) =2^7 ^ (< = ° ) = ~2^ (6) 

where the largest component among the two is the one along y axis, because of the asymmetry of v z0 with respect to 
y axis. Conversely, v z q is symmetric with respect to the x axis (see eqs. (JH),©) and its partial derivative with respect 
to y ought to vanish in y — for any x. 

It is worth pointing out that the evolution equation for the classical vorticity uj in the relativistic case is more 
complicated than in non-relativistic fluid mechanics. However, it is still true that a non-vanishing initial value of 
the classical vorticity makes the fluid motion a vorticous one in general. This stems from the Carter-Lichnerowicz 
equation of motion (equivalent to the Euler equations) for a perfect fluid with one conserved charge [ic| |: 

(d^hu u - d u hu„) = Td„s (7) 

where h — (p + p)/n and s = s/n are the enthalpy and the entropy densities normalized to the charge density n, p is 
the pressure and T is the temperature. The vorticity tensor is usually defined as: 

= (d^huy - djiu^) (8) 

and the vorticity vector as [10]: 

It is quite straightforward to show from eq. ^ that the vorticity vector field has the non-relativistic limit: 



-e^^fU = -e a ^ Ufl {d v u a - cU„) (9) 



0,ivxvj (1.0) 
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so that w is a proper relativistic generalization of the classical vorticity uj = (1/2) V x v. The time component of a; is: 

iv° = ^7 2 v • V x v 

Thus, if V x v / then u> ^ and f2 ^ 0. It can be shown, starting from ([7]), that the spacial part of the vorticity 
tensor, that is f2 = V x /i 7 v, fulfills the Helmholtz vorticity equation: 



^ = V x (v x n) 



(11) 



provided that the fluid is isentropic, i.e. Vs = 0. All classical consequences of the vorticity equations then hold in 
relativity provided that uj is replaced by fl. 

Let us now study more in detail the fluid equations of motion at the time t = 0. We will write the Euler equation, 
instead of the Carter-Lichnerowicz, for a perfect ultrarelativistic fluid, with equation of state p = p/3. Accordingly: 



(p + p)(u- d)u» = g^dvp - (u ■ dp)u» 
and focus on the transverse components at the time t = 0, when u x = u y = 0, namely: 
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for i = 1, 2 where the equation of state has been used. Multiplying both sides by 7q and manipulating the derivative 
on the right hand side: 
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Now, since 7 2 = 1/(1 — u 2 ) : 



and eq. (JT^J) becomes: 
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Two terms are then responsible for the initial transverse velocity increase: the first is related to the gradient of the 
energy density in the observer frame, while the second depends on the gradient of the initial velocity field, i.e. on 
the vorticity according to eq. ([6]). If v z q was independent of x,y the second term would vanish and the transverse 
expansion would then be driven by the energy density gradient only, like in the usual picture. Because of the eccentrity 
of the overlap region (see fig. [5]), the system gets an initial kick larger in the x direction than in y and an anisotropy 
in the final spectra ensues. If, however, the second term is included, the expansion gets an additional contribution 
because dv z o/dx < (see fig. |4]), and v z q is negative for x > and positive for x < so that altogether the second 
term drives an increase of u x for x > and a decrease (starting from zero) for x < 0. Moreover, the expansion rate 
related to this term will be larger in the x direction than in y because, expectedly, dv z o/dx > dv z o/dy, (see discussion 
following eq. ©), thereby enhancing the elliptic flow. In other words, besides the geometrical anisotropy, elliptic flow 
gets a finite contribution from an initial kinematical anisotropy of the longitudinal velocity. This enhancement of 
elliptic flow can be seen as a centrifugal effect owing to angular momentum conservation: particles with a momentum 
orthogonal to J, i.e. directed along the reaction plane get an additional momentum kick with respect to those emitted 
along J. 

It is now interesting to make an estimate of how large this contribution is in our simple scheme. Assuming that the 
energy density is proportional to the total energy of nucleons in the overlap region so that: 



(Po +Po)lo -Po = 



2 ' ' ^ =±-\T{x-b/2,y)+T{x + b/2,y)]^M» 



(17) 



Az dxdy Az 

and, by using ([5]) and the equation of state p = p/3 we can obtain the expressions of the initial proper energy density 
P- 
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FIG. 5: Initial longitudinal velocity profile along the reaction plane y = for two different impact parameters for the collision 
of two hard sphere nuclei with radius 7 fm. 



and the flow velocity v z q: 




which is shown in fig. [5] for the case of hard sphere nuclei with radius 7 fm. According to eq. (rho), the proper energy 
density is an even function of x, as it was expected with the assumption (|17p . while v z q is an odd function of x. Also, 
tt can be seen from fig. [5] that v z o has a singular derivative at the edge of the overlap region which is due to the hard 
sphere assumption; such singularities disappear with smooth density profiles. By using (fT9|) , (fT8|) , and (fT"7|) we 
can compute the ratio of the second to the first term in eq. (116|) for the x axis: 
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(20) 



in order to evaluate the importance of the vorticity term for the expansion rate. This ratio is shown in fig. [5] for the 
case of hard sphere nuclei for two different y values at an impact parameter 6 = 6 fm. It is seen that the second term 
is a consistent fraction of the first term even near the collisions centre x = (about 20%) while it steeply increases 
at larger x values; at the boundary of the x interval the ratio shows spikes owing to the sharp sphere assumption 
and it is not shown. Of course, these numbers refer to an oversimplified example and just for the initial expansion 
kick, but the conclusion that the longitudinal velocity gradient cannot be neglected in more realistic hydrodynamical 
calculations should hold. 

As has been mentioned, in some hydrodynamical calculations [H, [llj, a non- vanishing angular momentum of the 
plasma is tacitly introduced by enforcing an asymmetric x dependence for the proper energy density in peripheral 
collisions keeping the Bjorken longitudinal scaling, i.e. the independence of v z on the coordinates x,y. Thereby, 
the longitudinal momentum density ([5]) conservation is fulfilled even though v z is independent of x and the angular 
momentum conservation Q is also fulfilled. We think that this assumption is quite unnatural. Firstly, it cannot 
hold in our specific example of instantaneous thermalization at infinitely large energy (with infinitesimally thin fluid 
in fig. [J) because the only velocity which is compatible with symmetry and independent of x is 0, thus making both 
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momentum and angular momentum density vanishing. However, even in the more realistic and more general case 
of finite thcrmalization time, it does not lead to the same flow velocity field as in the case of angular momentum 
conserved through Bjorken scaling breaking because of the absence of the vorticity term. This can be shown by 
enforcing the equality of angular momentum densities in the two approaches: 



4- .2- 4 2 



(21) 



where quantities with a tilde on the left hand side are such that only p depends on x while on the right hand side we 
have the standard ones in our approach. From the above equation follows: 
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Using (f2"2"|) and (|2"Tj) to obtain dp/dx in the equation of motion at the time t — (|13[) . we get, after some manipulations: 
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Conversely, in an approach where velocity is uniform with the same angular momentum density, one would have: 

1 dp 



(23) 



du x 
~~dt 



t=o 



4po7o dx 



t=o 



(24) 



Therefore, even though the same angular momentum density was enforced by modifying the energy density profile, 
the expansion rate could be consistently different from the one with non- vanishing vorticity because of the factor 7/7 
and, chiefly, the additional term proportional to the derivative of longitudinal velocity which in general speeds up 
expansion, as we have seen. 

The ratio of the expansion-driving terms in eq. I|20p depends only on geometry and not on the centre-of-mass 
energy, because so do both p and v z0 according to eqs. (fT5| . (fH?)) and the expressions (|5|) . (fT7|) . This apparent energy- 
independence is just a specific feature of our simple scheme where the longitudinal dimension was shrunk to a very 
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FIG. 7: Elliptic flow coefficient v^ 7 ^ as a function of pr for hadrons originated from a spherical spinning plasma at a chemical 
freeze-out T=165 MeV and radius 10.1 fm for ui/T = 0.03. The elliptic flow would simply vanish if J = 0. 



thin slab. In fact, this cannot be the case in a more realistic one, where thermalization is not instantaneous and, 
therefore, the gradient of v z is distributed on a larger volume. In other words, the relative vorticity contribution will 
not be as large as it turned out to be by enforcing instantaneous thermalization in an infinitesimal slab Az. This can 
be better seen by rewriting the angular momentum, for a perfect fluid, as: 



x 2 fx 2 
d 3 x xx tv = I d 3 x Vy x /ry 2 v = / d 3 x V x y ^7 2 v - 



d 3 x — V/17 2 xv-/ d 3 x x 2 h^ 2 uj 



d 3 x — V x hrfv 
2 ' 



(25) 



where Wi = T = hj 2 Vi is the momentum density, h = p + p is the enthalpy density and u> — (1/2)V x v; in the above 
equation we assumed that the enthalpy density vanishes outside a compact region. In general, the sum of the two 
terms in eq. (|25|) is constrained by angular momentum conservation but their relative contribution to it can, and will, 
vary with the centre-of-mass energy. Most likely, at lower energy, the vorticity term will be less important whereas, 
at higher energy, its relative contribution should approach the limiting one calculated in our simple scheme because 
thermalization is expected to be faster, with an initial denser plasma and a higher angular momentum density. If this 
is the case, at the LHC, a further increase of the elliptic How with respect to RHIC ought to be observed. 



IV. ELLIPTIC FLOW FOR A SPINNING SYSTEM 



Taking viscosity into account implies a modification of the hydrodynamical equations but this should not affect 
our conclusions. This can be understood with a simple argument: angular momentum has to be conserved anyway 
and dissipative effects such as viscosity will speed up entropy increase. Thus, the system will tend to the maximal 
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entropy configuration which, for a system with finite angular momentum and finite volume, is a rigidly spinning fluid, 
with velocity field v = 10 x x, lo being a constant vector related to the total angular momentum J @, [13] • Of course, 
the quick expansion will prevent the system from reaching global equilibrium before decoupling, still the system will 
evolve towards that configuration. 

Hence, although viscosity, along with other dissipative forces, will not be able to create a fully equilibrated rigidly 
spinning plasma fireball, it is interesting to show that in this ideal case the effect of a finite angular momentum on 
elliptic flow and other observables is remarkable. 

That a globally spinning interaction region brings about an anisotropy in the particle azimuthal spectra was argued 
many years ago by Hagedorn [l3T ] and recently rediscussed in ref. [bj ]. Assuming statistical hadronization for a fully 
equilibrated subsystem of the plasma, the elliptic flow coefficient can be calculated as that of a rigidly spinning ideal 
hadron-resonance gas with angular momentum such that J u < J and fixed angular velocity u> = (1/2) V x v 
parallel to it and linked to J w by a thermodynamic relation which is linear for small lo/T values 0. In the Boltzmann 
approximation for primary hadrons, this reads Q: 



(26) 



(j) _ v /l " |a ' xx| 



Ki(m T i /l-|wxx|? IT) 



where Ki, I„ are modified Bessel functions and T is the global temperature. It should be stressed that the global 
temperature T in a spinning relativistic gas is related to the local proper temperature To, measured by a comoving 
thermometer, by the relation Q: 

T °( r ) = h T 5 a W 

VI — 10 T 

where r is the distance from the rotation axis u). Since it is the local, and not the global, temperature which determines 
the phase of the system, the decoupling should occur when the highest local temperature reaches the critical value T c 
for the quark-hadron transition, that is when: 

- T = T c (28) 
y/1 - oj 2 R 2 V ' 

being R the maximal distance from the rotation axis. 

The behaviour of the "rotational" (which would simply vanish if J = 0) as a function of px for primary hadrons 

is very similar to that driven by pressure gradients in usual hydrodynamical calculations, and it turns out to be almost 

independent of the particle mass. It is shown in fig.[7]for lo/T = 0.03 at the chemical freeze-out temperature T c = 165 

MeV, for a spherical source with radius R = 10.1 fm, for a total angular momentum J w ~ 10 , i.e. of the same order 

of the J of the interaction region at RHIC energies. The v ( 2 J) of primary hadrons from a globally spinning region 

would be therefore very large, although resonance decays should lower the final one consistently. 

V. POLARIZATION 

Elliptic flow is not a unique consequence of an intrinsic rotation. There is, however, a distinctive signature thereof: 
a polarization of the emitted hadrons along the angular momentum direction (in the observer frame) . That a large 
angular momentum in peripheral heavy ion collisions could give rise to polarization of the final hadrons has been first 
proposed in refs. @, where a quantitative assessment was performed within a perturbative QCD framework, with the 
polarization of quarks assumed to be effectively transferred to final hadrons. Recently, it has been pointed out that a 
plasma with polarized quarks could be probed by observing the polarization of direct photons [l5| . We take a different 
approach here and we determine the polarization of particles invoking local thermodynamical equilibrium and the 
statistical hadronization dogma which is succesfull in describing hadronic multiplicities: every multihadronic state 
compatible with conservation laws is equally likely. Therefore, since the total angular momentum is not vanishing, 
when the plasma hadronizes, available spin states will not be evenly populated and a net polarization of the produced 
hadrons will show up. In this approach, there is no need to invoke any special dynamical mechanism for the polarization 
of quarks to be transferred to hadrons, as it should happen as a consequence of the statistical nature of this process. 

The proper polarization vector IIq of particles in a relativistic rotating ideal gas has been calculated by the authors 



^ 1 , w 
Iln = — tanh — 
2 2T 



e „ lj pp 

LO — 



m m(e + m) 



(29) 
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for spin 1/2 particles and: 



n, 



L^dn—- 



, ne 



j/T 



— W 

m m(e + to) 



(30) 



for generic spin particles, where e is the energy and p the momentum of the particle. The polarization along o>, 
i.e. IIo • Cj turns out to be maximal for particles emitted orthogonally to u) (i.e. along the reaction plane for an 
equilibrated spinning system) and increases for increasing pr up to momenta of the order of 2mT/uj where the 
rotational grand-canonical ensemble scheme fails and more complicated expression arise Q. Also, the vector mesons 
show spin alignment in that the 00 component of the spin density matrix turns out to be different from 1/3 and 
specifically [9(: 



which, for small co/T, reduces to: 



2cosh(cj/T) + 1 



cosh(w/T) 



p 2 w 2 



(1 -cosh(w/T)) 



Pujoo(p) 



1 l-3(p-w) 5 



18 



J^2 



(31) 



(32) 



It is interesting to note that the polarization (more generally the spin density matrix) depends on the ratio between 
angular velocity and global temperature, that is, using ([27]) on juj/T , being T the local temperature. Therefore, 
it can be conjectured, by invoking locality principle, that a polarization should appear in a generic accelerated 
hydrodynamical cell at local equilibrium fully determined by local quantities. Hence,u; is to be plausibly replaced by 
the vector: 



1 



(33) 



which is the local angular velocity for a general trajectory, according to the Frenet formulae. If this conjecture is true, 
every hadronizing hydrodynamical cell will produce hadrons with polarization vector (|30p with u> equal to the right 
hand side of (|33|) . 

The expected polarization values are of the order of cj/T, which is reasonbly some percent or less (see previous 
section) but they should increase with particle momenta up to momenta of few GeV's and hopefully become observable 
1 . However, the expressions (I3TH) and (|3"Tj) refer to primary hadrons, i.e. those emitted from the source at decoupling 
and resonance decays can further dilute the polarization, so that a more detailed study is needed. 

It is difficult to predict the evolution of polarization values as a function of centre-of-mass energy. However, it 
can be argued that they should increase by considering that angular momentum density at freeze-out increases as a 
function of •s/s/v./v- This happens because the size of the system at freeze-out increases approximately logarithmically 
whereas the angular momentum of the interaction region increases linearly with ^/s NN (see eq. (|3])). Since the angular 
momentum density must be somehow related to the final local angular velocity (|33|). a fair conclusion follows that 
polarization effects should increase with the collision energy. 



VI. SUMMARY AND CONCLUSIONS 



We have pointed out that angular momentum conservation in peripheral ultrarelativistic heavy ion collisions at 
very high energy should give an additional contribution to the azimuthal momentum anisotropy, thereby enhancing 
the elliptic flow coefficient v^- By using a very simple hydrodynamical scheme, we have shown that taking angular 
momentum conservation properly into account implies, most likely, a non-uniform longitudinal flow velocity in the 
transverse plane breaking the usual assumption of Bjorken scaling. This in turn generates a non-vanishing initial 
vorticity term in the equations of motion which enhances the transverse expansion rate and may be able to balance 
the elliptic flow deficit observed by Heinz and Song [j| in minimally viscous hydrodynamical calculations. Angular 
momentum conservation is also implemented in current hydrodynamical calculations by keeping the Bjorken scaling 
hypothesis, but the resulting expansion rates are different. We expect this effect to be more visible at very large 



1 Recent measurements by RHIC experiments set a limit on average A polarization to be ~ 0.02 [l6| 
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energies where the vorticity contribution to the angular momentum density tends to an upper geometrical limit, that 
we have analyzed in our simplified scheme. Hence, we predict that i> 2 should increase from RHIC to LHC energy, 
although we cannot give a definite quantitative estimate. 

The most characteristic signature of the vorticity induced by angular momentum conservation would be a polariza- 
tion of the emitted particles, which is predicted to be, in the observer frame, for a globally spinning system, orthogonal 
to the reaction plane and maximal for particles with momentum parallel to the reaction plane in case of a globally 
spinning plasma. A quantitative assessment of these effects for the actual hydrodyamical evolution is very difficult, 
but we argued that the polarization should be there for a general accelerated fluid motion. Also this effect should be 
better observed at the LHC, where the angular momentum density should be larger. 
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